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Abstract 



We study the factorization scheme dependence of the next-to-leading order 
inclusive one jet cross section da /oIEt- The scheme is varied parametrically 
along the direction that transforms the MS scheme to the DIS scheme: we 
introduce a parameter A such that A = is MS and A = 1 is DIS. The 
factorization scale \x is also varied. We observe a change of ±9% in the cross 
section for Et = 50 GeV when \i and A are varied in the range — 2 < A < 2, 
E T /8< fi< 2E T . This grows to ±32% for E T = 400 GeV. 
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Calculations in perturbative QCD involve conventions. For example, one conventionally 
uses the MS scheme for renormalization. Within the MS scheme, one needs to specify 
a renormalization scale //uv- If one does a calculation at next-to-next-to- leading order 
(NNLO), it is possible to vary the renormalization scheme by setting the third coefficient 
in the (3 function, 6 3 . It is useful to study the dependence of the calculated cross sections 
on the chosen conventions. Thus, with a NNLO calculation of the total cross section for 
e + + e~ — > hadrons, one can examine the dependence of the cross section on /iuv and 63 [|TJ. 
Studies like this have a dual purpose. First, we learn about the extent to which our arbitrary 
choices affect the results. Second, we obtain an estimate, or at least an approximate lower 
bound, on the theoretical uncertainty in the calculated order cross section arising from 
uncalculated terms of order a^ +1 and higher. This is because among the uncalculated terms 
are contributions containing factors of In (/iuv) and 63 that cancel the /iuv and 63 dependence 
in the calculated terms. We cannot find all of the a^ +l terms (since we lack the calculation), 
but at least we know the size of some of them. 

In cross sections involving initial state hadrons, the theoretical formula contains parton 
distribution functions as factors. The calculated result at order depends on the conven- 
tion we use for collinear factorization. Put another way, the result depends on the definition 
of the parton distribution functions. In particular, the result depends on the factorization 
scale /x co n (the scale that appears in the parton distributions f a /A{ x i AWi))- I n NLO cal- 
culations of hard cross sections with initial state hadrons, it is common to examine the 
dependence of the result on /x co ii - 

In this paper, we examine a bigger space of convention dependence. The dependence of 
the parton distribution functions on the factorization scale has the form 

fa/A(x,P) = fa/A{x,fl) 

+ 0{a% (1) 

Here K ® b (£) is the lowest order Altarelli-Parisi kernel. If we like, we can redefine the parton 
distributions according to 

fa/A&Ll) = fa/A(,X,fi) 

+ o^/i) \j f 1 ^ ^tf(^)A /A (^). (2) 
Ztt j=q Jx t, b 

Here K^(^) is still the Altarelli-Parisi kernel, so that if we put Ao = ln(/i 2 //! 2 ) then the 
J = term amounts to a change of scale, considered at lowest order in a s . The essential 
new ingredient is that we now have added more terms with kernels (£), which can be 
anything we please. 

Let us suppose that we are calculating a cross section a. If we calculate the hard 
scattering (partonic) cross section at leading order, then the order a s change in definition 
of the parton distributions caused by varying one of the A j will lead to a relative order a s 
change in the cross section, 
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-^-(xaxa s (LO). (3) 

If we calculate the hard scattering cross section at N fc LO, then the N fc LO contribution to the 
hard scattering functions will contain contributions proportional to the Aj raised to powers 
up to k, so that the Aj dependence is partially canceled and 

da 



— — oc a x a. 



k+l 



(N fc LO). (4) 



Suppose that we have a calculation at NLO = N 1 LO. Then the change in the cross section 
induced by changing some of the Aj provides an indication of uncalculated contributions 
to a that are suppressed by a factor a 2 . Looked at another way, changing some of the 
Aj provides an indication of the sensitivity of the calculation to the arbitrary choice of 
factorization convention. 

Note that the question of the dependence on the factorization convention is analogous 
to the question of the dependence on the renormalization convention. However, in the case 
of renormalization we are really varying one quantity, a s , so that there is a one parameter 
space of conventions to explore if we work at first order, Aa s oc a s x a s . In the case of 
factorization, we are varying functions, the parton distribution functions, so that there are 
an infinite number of parameters to examine at first order, Af oc / x a s . 

In this paper, we restrict the parameter space by taking N = 1, with being the 
kernel such that Ai = 1 corresponds to changing from the MS scheme to the DIS scheme. 
We examine the cross section da/dE? for pp — ► jet + X calculated at NLO and look at how 
this cross section depends on A and Ai. 

In our calculation, we eliminate the J = term in Eq. (|2|) and instead change the scale 
/i in the parton distributions directly, using the parton evolution that is incorporated in the 
CTEQ4M set of parton distributions.^] This differs from using the J = term at fixed \i 
by terms of order a 2 . Thus, we have a two dimensional parameter space of conventions to 
explore, specified by \i = /x co n and A = Ai- 

Using the program 0, we computed the one jet inclusive cross section da/dEx, where 
Et is the transverse energy of the jet. Suitable modifications for implementing changes in 
scheme were made. Jet definitions followed the Snowmass cone algorithm with the cone 
radius R = 0.7. We averaged over the rapidity y in the range 0.1 < \y\ < 0.7. The 
renormalization scale was fixed at /xyv = -Et/2 to limit the number of variables. The 
factorization scale fi was varied in the range Et/8 < fi < 2Ex- The factorization scheme 
parameter A was varied in the range — 2 < A < 2, with A = corresponding to the MS 
scheme and A = 1 to the DIS scheme. 

We calculate the DIS parton distribution functions from the MS ones, with CTEQ4M as 
our standard starting distribution: 

Ux,^X) = f^(x^) + X^ £i^W£)/P(^)- (5) 

27T Jx 4 , 



1 The CTEQ4M parton distributions incorporate evolution with both a s and a 2 s terms in the 
evolution kernel. 
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The kernel consists of the usual DIS functions 0, 

KM,) = -K!}n,) = c,. 1,1(1 "-'^ 3 



2n f ~ 99 



KV}(x)=Ti 




[1 + x) ln(l — x) 
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-8x 2 + 8x-l 



(6) 



with q = u,u,d,d, ... Here, Cf = 4/3, Tp = 1/2 and the other are equal to zero. 
The jet cross section is originally defined as 
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^] f dx A f dx B y\f a {x A )f b {x B ) 
2nJ Jo Jo 



a. 



a (a, b, x A , x B ) + — a x (a, b, x A , x B ) 

Z7T 



(7) 



where &o is the leading-order (or Born) term, and o\ is the next-to-leading-order term. We 
can selectively turn on or off the NLO term in the cross section. We reexpress this cross 
section in terms of the new parton distributions f(x), dropping contributions suppressed by 
two powers of a s compared to the Born cross section, 



a 



^) dx A dx B ^2\f a (x A )f b (x B ) 



a.b 



& Q (a,b, x A ,x B ) + — di(a,&, x A ,x B ) 
In 



xa (a,b, x A ,x B ) \, 



d£ 



(8) 



where Ao = ln(4/x 2 /i?|,) and ot s = a s (/xuv)- The added terms modify &i from the standard 
MS scheme, that is the MS scheme where \l = Et/2. (Of course the J = term is already 
included in the standard program 0.) Changing both parton distribution and jet cross 
section definitions in this way exactly cancels out contributions of order a s to the change in 
a. What remains are the order ct\ contributions and the extent of the change in the cross 
section gives us an estimate of these NNLO terms. 

We first calculate the jet cross section da/dEx at Born level at the transverse energy 
Et = 100 GeV. Fig. |I|(a) shows a contour plot of the fractional difference of this computed 
cross section from the standard MS cross section when A and fi are varied. The contours 
represent 1% changes in the cross section. The Born level cross section is very sensitive to 
changes in conventions; it varies by 40% from its standard value in the region shown in the 
plot. 

In Fig. |l](b), we show the fi and A dependence of the cross section calculated at NLO. The 
NLO terms cancel out most of the convention dependence and a saddle region is observed. 
The maximal change from the standard cross section is only 8%. There is a rather broad 
region in the parameter space under consideration where the NLO cross section changes 
little. 
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In Fig. 0, we investigate how this result depends on Et in the range 50 GeV < Et < 
400 GeV. The sensitivities of both the Born cross section (not shown) and the NLO cross 
section to \x and A increase with increasing Et- However, the Born cross section is always 
far more sensitive to the A and \x parameters than the NLO cross sections. At NLO, saddle 
regions are observed for all values of Et shown, although the saddle point is not always at 
the same position in the A /i plane. 

We can offer two comments about the results shown in Figs. [I] and 0. First, the choice 
/i = Et/2 is quite felicitous: with this choice, the dependence on A is almost nil. Second, 
there is a lot more convention dependence at Et = 400 GeV than at Et = 100 GeV or 
Et = 50 GeV. (One misses seeing this if one looks only at the \i dependence at A = 0. 
Instead, one needs to vary A and \i simultaneously.) 

One can use the results of Figs. [I] and |2] to estimate a theoretical error to be ascribed to 
the calculation, as described in the introductory paragraphs. We can estimate the error as 
the amount by which the cross section changes when one makes a "substantial" change in /i 
and A. But we need an a priori decision about what range of \i and A represents a substantial 
change. Evidently such a decision must be subjective but we cannot avoid making a choice. 

Consider /i first. This parameter is supposed to represent the typical momentum that 
flows in loops of Feynman graphs for the process. For instance, a good guess would be 
H = Et sin(i?) where R = 0.7 is the cone size. Thus \x = Et/2 is not an unreasonable 
guess. Clearly such an estimate cannot be good to much better than a factor of two, but our 
intuition is that the estimate should not be off by much more than a factor of two either. 
Thus we adopt a factor of two as a measure of a "substantial" scale change. 

Now consider the parameter A. The choice A = represents the MS convention. The 
choice A = 1 represents the DIS convention, in which the contribution to deeply inelastic 
scattering from gluon initial partons is canceled by choice of convention. Since this represents 
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a qualitative alteration in how deeply inelastic scattering appears to happen, it seems to us 
that AA = 1 represents a "substantial" change. 

We thus propose that the variation in the cross section when \i and A vary in the range 
{Alog 2 (/i) = ±1,AA = ±1} represents a minimum theoretical error in the sense that it 
would be a surprise (to us anyway) if the difference between the NLO result and a NNLO 
result were less than that. A more conservative error estimate would come from doubling 
the range: {Alog 2 (/i) = ±2, AA = ±2}. 

With this understanding, we can read error estimates off of Fig. |^. The minimum error 
estimate varies from 3% at = 50 GeV to 7% at = 400 GeV. The conservative error 
estimate varies from 9% at E T = 50 GeV to 32% at E T = 400 GeV. 

Klasen and Kramer have carried out an interesting investigation Q that is in some 
respects similar to that reported here. They calculated the jet cross section at large Et 
using the CTEQ3D parton distributions, which are defined with the DIS convention. They 
compared this cross section to the cross section calculated using the CTEQ3M parton distri- 
butions, which are defined with the MS convention. In each case, they used the a appropriate 
to the parton convention. They found a large difference. We agree with this result. We 
have calculated da/dE? at E? = 400 GeV with fi = Et/2 using CTEQ4D parton distri- 
butions and the DIS version of a. We find that this DIS cross section is 33% greater than 
the corresponding MS cross section calculated using CTEQ4M parton distributions and the 
MS version of a. How can this be consistent with Fig. [|, which indicates that the DIS cross 
section (A = 1) differs from the MS cross section (A = 0) by less than 1% for \x = Et/27 
The answer is that the DIS parton distributions obtained from the CTEQ4M distributions 
by using the transformation (|5|) are not the same as the CTEQ4D parton distributions. 

We can comment in more detail. The CTEQ4D distributions are not obtained from the 
CTEQ4M distributions by using Eq. @ but are instead obtained by fitting the world's data 
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FIG. 3. CTEQ4M, CTEQ4D and calculated gluon distribution at // = 200 GeV 

using a DIS version of NLO theoretical formulas. In Fig. |3] we display the gluon distribution 
at /i = 200 GeV as given by the CTEQ4M set, by the CTEQ4D set, and by the calculated 
transformation (|5|) from the CTEQ4M set. We see that the DIS gluon distribution calculated 
using Eq. ([5]) is negative for x > 0.5 for fi = 200 GeV. The CTEQ4D gluon distribution is 
positive for all x. (This is a constraint imposed in the fitting procedure.) Thus the CTEQ4D 
gluon distribution is substantially larger at large x than it might have been. This does not 
create a bad fit since there is essentially no data that constrains the gluon distribution at 
large x 0. The larger gluon distribution leads to a larger jet cross section. Thus the 33% 
change in the calculated cross section can be attributed to the uncertainties in fitting parton 
distributions, arising ultimately from the fact that the gluon distribution at large x is not 
constrained by data. 

In summary, we have investigated a two dimensional space of factorization schemes. The 
dependence of the jet cross section on the two parameters [i and A is moderate in the range 
50 GeV < Et < 400 GeV. The sensitivity to fi and A increases as Et increases. Within 
the space investigated, the choice fi = Et/2, A = (the MS scheme with a standard choice 
of scale) seems as good as other nearby choices. It will be interesting to see whether these 
conclusions change when we investigate a bigger space of factorization schemes. 
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